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1. Introduction 



The 'running' of the renormalised quark mass as the scale M is changed is controlled by the j8 
and 7 functions in the renormalisation group equation, defined by 



j3V(M)) 



dg s (M) 



dlogM 



dlogM 



(1.1) 
(1.2) 



where the bare parameters are held constant. These functions are given perturbatively as power 
series expansions in the coupling constant, 



(1.3) 



j8 S (g) = -b g 3 - byg 5 - b^ 1 - b\g 9 

£ {g) = d m og 2 + d s ml g A + df n2 g 6 + d s m3 g s + 



The first two coefficients of the j8 -function and first coefficient of the y m function are scheme 
independent, 

1 / 2 \ 1 



ba 



(47T) 2 i n 3 



(4n) 4 



38 

102- —n f 



(1.4) 



and 



{An) 2 ' 



(1.5) 



while all others depend on the scheme chosen. 



We may immediately integrate eq. ( 1 . 1 ) to obtain 



M 



— s = [b g s {M) 2 Y%zx V 



2b gS(MY 



exp 



8 s (M) 







+ 



m) bot 3 bit; 



(1.6) 



The renormalisation group invariant (RGI) quark mass is defined from the renormalised quark 
mass as 



mf = AZi(M)m s (M) = AZ s m {M)Z s m (M)m b q ar ' = Z™m b q a " , (1.7) 



where 



[AZ^(M)]- 1 = [2b g S (M) 2 ] 2 ™o exp . 



S(M) 



TO) M 



(1.8) 



and so the integration constant upon integrating eq. ( |j~T| ) is given by A s , and similarly from eq. ( 1 .2) 
the integration constant is m^'. A s and m q u are thus independent of the scale. (Note that although 
the functional form of AZf n (M) is fixed, the absolute value is not; conventions vary for its defini- 
tion.) Also for a scheme change S — > S' (it is now sufficient to take them at the same scale) given 
by 

/ = G(/)=/(l + I ?1 (/) 2 + ...), (1.9) 



Analogous definitions hold for other quantities which depend on the scheme and scale chosen. 
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rrvV remains invariant, while A changes as A s ' = A s exp(?i / (2bo)). Note also that analytic expres- 



sions for the integrals in eq. ( |L7[ ) or eq. ( |1.8| ) can be found for low orders, for example to two loops 
we have 



AZ? n (M) = [2b (g s (M)) 



21 255 



1 + !V(M)) 2 
bo 



b n d S , - 
U ml 



(1.10) 



Thus we have a convenient splitting of the problem into two parts: a number, m", which involves 
a non-perturbative computation, and is the goal of this paper and, if desired, an evaluation of AZ^ 
which allows the running quark mass to be given in a renormalisation scheme S. 



2. Simulation 

We have estimated the light quark masses in the MS scheme at 2GeV by first using the axial 
Ward identity (AWI) to determine the lattice quark mass. This is renormalised using the RI' — MOM 
scheme [|I|] (for our variation on the method see |^]), converted to a RGI form as described in 
section 1 and after the continuum limit has been taken rewritten in the MS scheme. Further details 
and results are given in [Q]. We perform our simulations with two flavours of non-perturbatively 
clover-improved dynamical Wilson fermions and Wilson glue. Using these actions, the QCDSF 
and UKQCD collaborations have generated gauge field configurations with the parameters given 
in Table [l]. We also use configurations generated by the DIK collaboration which have been made 
available through the ILDG. This large set of lattices enables us to extrapolate to the chiral and the 
continuum limit. 



3. Quark masses 

As a first check we perform the chiral extrapolation for pseudoscalar mass. In Fig. [I] we 
plot (am ps ) 2 against am^ WI together with the fit result. Our data shows that {am ps ) 2 goes to at 
=0. This means that even at our lightest quark mass the data is not seriously effected by 
either an Aoki phase or weak 1st order phase transition. 




AWI 

am 



Figure 1: (am ps ) 2 versus am™ 1 for j3=5.29. 
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Table 1: Overview of our lattice parameters. For the translation into physical units the Sommer 
parameter^ with ro = 0.467 fm (see [Q] and Q|) has been used. 





K 


N 3 xT 


m % [GeV] 


a [fm] 


L [fm] 


Ntraj 


5.20 


0.13420 


16 3 x 32 


1.007(2) 


0.115 


1.8 


^(5000) 




0.13500 


16 3 x 32 


0.833(3) 


0.098 


1.6 


^(8000) 




0.13550 


16 3 x 32 


0.619(3) 


0.093 


1.5 


^(8000) 


5.25 


0.13460 


16 3 x 32 


0.987(2) 


0.099 


1.6 


^(6000) 




0.13520 


16 3 x 32 


0.829(3) 


0.091 


1.5 


^(8000) 




0.13575 


24 3 x 48 


0.597(1) 


0.084 


2.0 


^(6000) 


5.29 


0.13400 


16 3 x 32 


1.173(2) 


0.097 


1.6 


^(4000) 




0.13500 


16 3 x 32 


0.929(2) 


0.089 


1.4 


^(5600) 




0.13550 


24 3 x 48 


0.769(2) 


0.084 


2.0 


^(2000) 




0.13590 


24 3 x 48 


0.591(2) 


0.080 


1.9 


^(5000) 




0.13620 


24 3 x 48 


0.395(3) 


0.077 


1.9 


^(3000) 




0.13632 


32 3 x 64 


0.337(3) 


0.077 


2.5 


47(1000) 


5.40 


0.13500 


24 3 x 48 


1.037(1) 


0.077 


1.8 


^(4000) 




0.13560 


24 3 x 48 


0.842(2) 


0.073 


1.8 


^(3000) 




0.13610 


24 3 x 48 


0.626(2) 


0.070 


1.7 


^(4000) 




0.13640 


24 3 x 48 


0.432(3) 


0.068 


1.6 


^(3000) 



We use the next to leading order (NLO) chiral perturbation theory (#PT) to estimate the quark 
masses, 

romf = cf [(r m K+ ) 2 + (r m^) 2 - (r m^) 2 ] 

+ (cf - cf) [(r m K+ ) 2 + (r m K0 ) 2 ] {r m n+ ) 2 
+i(cf + cY)[(r m K+ ) 2 + (rom K ») 2 ] 2 
-{cf + cf){r Q m n+ f 

-cf [(r m K+ ) 2 + (rom^o) 2 ] [(r m K+ ) 2 + (r m K o) 2 - (r G m n+ ) 2 ] 

x In ((r m K +) 2 + (r m K o) 2 - (r m^+) 2 ) 
+cf (r m„ + ) A \n(r m„ + ) 2 (3.1) 

rQm^ = cf(r m 7t+ ) 2 + (cf + cf)(r Q m !t+ ) 4 -cf( r Qm J[+ ) 4 \n(r m J[+ ) 2 . (3.2) 
The fit function to determine cf and cf 3 , i = b,c,d is 
r n m RGI 

, \ 2 = c a'+ CT (rom%) 2 + cf (romps) 2 + cf ' {(r m s ps ) 2 -2(r m ps ) 2 )\n(r m ps ) 2 , (3.3) 

where m ps , m ps are the valence and sea pseudoscalar masses respectively (both using mass degen- 
erate quarks, since we found the relevant quantities am ps and am q to differ by < 1% between the 
degenerate quarks case and the non-degenerate quarks case). The first term is the leading order, 
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LO, result in #PT while the remaining terms come from the next non-leading order, NLO, in #PT. 
We note that to NLO, we can determine cf and cf , i = b,c,d using mass degenerate quarks and 



then simply substitute them in eqs. < p.l[ |3.2| ) 



To reduce the total error on the result, it proved advantageous to use eq. ( |3.1| ) to eliminate 



from eq. (3.3) in terms of 



c„ 



r m" 



(r m K+ ) 2 + (r Q m K o) 2 - (r w%+) 2 



(3.4) 



This results in a modified fit function of the form 

-K3 



{romps 



c 



«3 



+ cf [{rQtn ps ) ~ d b ) + cf \{r m ps ) ~ d c 



+ cf [{(r m s ps ) 2 - 2(r m ps ) 2 ) ln(r m p ,.) 2 - d d ] 



(3.5) 



where di (i = b,c,d) can be read-off from eq. ( |3.1[ ) and have the effect of shifting the various terms 
in the fit function by a constant. 

0.16m ■ 1 ■ 1 ■ 1 ■ 



0.14- 
0.12- 
0.10- 
0.08- 
0.06- 
0.04- 
0.02- 

o.oo- 

-0.02 - 




0.01 



0.02 

w 2 



0.03 



0.04 



Figure 2: cf , cf 3 (z = b,c,d) and cf versus (a/ro) 2 . The open symbols represent the values of cf, cf 
(i = b,c,d) and cf in the continuum limit. 

In Fig. U we plot cf, cf (i = b,c,d) and cf against (a/ro) 2 . The coefficients of NLO are 
small compared to the LO coefficient. Finally, we find 



mf (2GeV) 



121(2)(3)(6)MeV for r = 0.5fm 
115(2)(3)(6)MeV for r = 0.467fm 



(3.6) 



where the first error is statistical and the second is systematic 3MeV. We have determined it from 
the effect on cf by changing the fit interval (rotn ps ) 2 < 5 to (rom ps ) 2 <4 or 6 or °°, i.e. include all 
the data. Furthermore the additional third (systematic) error is due to the uncertainty with which 
value to identify ro. 

For the light quark mass, we find that corrections from LO to NLO #PT are negligibly small. 
We shall just quote the LO result of 



<j(2GeV) 



4.57(05)(07)(23)MeV for r = 0.5fm 
4.34(05) (07) (23) MeV for r = 0.467fm 



(3.7) 
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Table 2: Comparison of the updated and previous values of m" s (2GeV), m"J(2GeV) and 
mf{2GeV)/m^(2GeV). 

previous new 
mf (2GeV) lll(6)(4)(6)MeV 115(2)(3)(6)MeV 

mf d {2GeV) 4.08(23)(19)(23)MeV 4.34(05) (07) (23 )MeV 

mf (2GeV)/m^(2GeV) 27.2(3.2) 26.6(1.8) 



where again the second and third errors are systematic. Finally, we see that the ratio 

mf (2GeV) 



<J(2GeV) 

4. Conclusion 



26.6(1.8). (3.8) 



We have updated our estimate for the light quark masses using results at smaller lattice spacing 
and smaller quark masses data. In Table ||, we compare the updated and previously published 
values [Hi. Our results are in rough agreement with other group's results. In order to improve the 
precision and accuracy of analysis, simulations of smaller quark masses and lattice spacing and 
2+ 1 flavours are needed. 
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